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Abstract 

Development of robust yeast strains that can efficiently ferment lignocellulose-based feedstocks is one of the 
requirements for achieving economically feasible bioethanol production processes. With this goal, several genes 
have been identified as promising candidates to confer improved tolerance to S. cerevisiae. In most of the cases, 
however, the evaluation of the genetic modification was performed only in laboratory strains, that is, in strains 
that are known to be quite sensitive to various types of stresses. In the present study, we evaluated the effects 
of overexpressing genes encoding the transcription factor {YAP1) and the mitochondrial NADH-cytochrome b5 
reductase {MCR1), either alone or in combination, in an already robust and xylose-consuming industrial strain of 
5. cerevisiae and evaluated the effect during the fermentation of undiluted and undetoxified spruce hydrolysate. 
Overexpression of either gene resulted in faster hexose catabolism, but no cumulative effect was observed with 
the simultaneous overexpression. The improved phenotype of MCR1 overexpression appeared to be related, at 
least in part, to a faster furaldehyde reduction capacity, indicating that this reductase may have a wider substrate 
range than previously reported. Unexpectedly a decreased xylose fermentation rate was also observed in YAP1 
overexpressing strains and possible reasons behind this phenotype are discussed. 
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Introduction 

Production of second-generation bioethanol from lignocel- 
lulosic biomass requires robust Saccharomyces cerevisiae 
strains with improved capacity to cope with the toxic com- 
pounds formed during the biomass pre-treatment, among 
which are 5-hydroxymethylfurfural (HMF), furfural, weak 
organic acids and phenolic compounds (Parawira and 
Tekere 2011). This has led to extensive studies to decipher 
mechanisms behind the compounds toxicity and the yeast 
natural tolerance responses to them and, among others, 
genes involved in detoxification and yeast tolerance to 
individual inhibitors have been identified, such as ADH6, 
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HAA1 or PMA1 (Haitani et al. 2012; Mira et al. 2010; 
Petersson et al. 2006); for a more exhaustive review, see 
(Almeida et al. 2009a; Liu 2011). YAP1 is another interest- 
ing candidate for industrial strain engineering because it 
encodes a transcription factor (Yaplp) that simultaneously 
controls a wide range of stress-related targets (Toone and 
Jones 1999). Notably, its overexpression has a beneficial 
role in the response of laboratory S. cerevisiae towards 
HMF, furfural, and different concentrations of hydrolysate 
(Alriksson et al. 2010; Kim and Hahn 2013; Ma and Liu 
2010; Sundstrom et al. 2010). Another interesting and 
complementary candidate for gene overexpression is MCR1 
that encodes the mitochondrial NADH-cytochrome b5 
reductase (Hahne et al. 1994; Meineke et al. 2008). Previous 
experiments performed in our group revealed that overex- 
pression of MCR1 in S. cerevisiae resulted in a reduced lag 
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phase and faster growth rate when the yeast was grown 
with high concentrations of acetic acid (Signori et al., 
personal communication). This weak acid is one of the 
inhibitors directly affecting xylose metabolism in S. cerevi- 
siae (Almeida et al. 2011; Bellissimi et al. 2009; Casey et al. 
2010; Helle et al. 2003). Still, considering that many of the 
studies about strain improvement towards hydrolysate-de- 
rived inhibitors concern laboratory strains, it is difficult to 
predict the real effect of these changes in an industrial, 
and more robust, strain background. 

The objective of the present study was to evaluate the 
effect of overexpressing YAP1 and MCR1, either alone 
or in combination, in process-like conditions, that is 
using a robust industrial S. cerevisiae strain and undetoxi- 
fied lignocellulosic hydrolysate. For this, the strain GSE16 
was chosen as background strain for engineering since it 
combines a robust industrial background with the ability 
to ferment xylose, using the xylose isomerase pathway 
(Demeke et al. 2013b). This robust strain was developed 
by a combination of different strategies including rational 
metabolic engineering, mutagenesis, evolutionary engin- 
eering, genome shuffling and meiotic recombination 
(Demeke et al. 2013a,b). Moreover, and as part of the 
genetic engineering strategy used in the current study, the 
overexpression of YAP1 was combined with the deletion 
of the chaperone-encoding gene APJ1 since deletion of 
this gene has been previously reported to enable growth 
on rich medium at inhibitory ethanol concentrations 
for the parental strain (Swinnen et al. 2012). The relevance 
of overexpressing YAP1 and MCR1 in an industrial 
background was confirmed for undetoxified spruce 
hydrolysate fermentation. We also uncovered unexpected 
interactions between YAP1 overexpression and xylose 
metabolism. 

Materials and methods 

Strains 

S. cerevisiae strains utilized in this study are presented in 
Table 1. Escherichia coli DH5a and E. coli NEB 5-alpha 
were used for sub-cloning and were grown on Luria- 



Bertani (LB) medium supplemented with 100 mg.L~ 
ampicillin, when required. Plasmids utilized in the study 
are described in Table 2. 

Molecular biology methods 

Standard molecular biology methods were used for all clon- 
ing procedures. (Sambrook and Russel 2001). Thermo- 
Scientific GeneJET plasmid miniprep kit (ThermoScientific, 
Lithuania) was used for plasmid extraction. E.Z.N.A Cycle- 
Pure Kit (Omega Biotek, USA) was used for purification 
of polymerase chain reaction (PCR) products. Qiagen 
Qiaquick gel extraction kit (Qiagen GmbH, Germany) 
was used to extract DNA from agarose gels. All DNA- 
modifying enzymes were purchased from ThermoScienti- 
fic. Primers for PCR and sequencing of DNA constructs 
were ordered from and performed by MWG (MWG-Bio- 
tech AG, Germany). All primers are shown in Additional 
file 1: Table SI. Transformation of E. coli was performed 
using the Inoue Method (Sambrook and Russel 2001). 
Either the lithium acetate (LiAc) method (Gietz et al. 
1995), a modified version of it that uses dimethyl sulfoxide 
(DMSO) (Hill et al. 1991) or electroporation (Benatuil 
et al. 2010), were used as transformation methods of 
S. cerevisiae. 

Construction of plasmids pJET1,2-attB-KanMX-attP, 
pJET1,2-attB-hph-attP and pBEVY-Nat-PhiC31 

The antibiotic markers KanMX and hph expressed 
under the TEF promoter and terminator were ampli- 
fied with the primers Fw-Al-attP-tefpr and Rv-A2- 
attB-teft and cloned into pJETl,2 (Thermo Scientific, 
Belgium) following the protocol of the kit. The result- 
ing plasmids were called pJETl,2-attB-KanMX-attP 
and pJETl,2-attB-hph-attP respectively. The PhiC31 inte- 
grase was amplified from pCMVInt (Addgene, Cambridge, 
Massachusetts, USA) using the primers Fw-PstI-PhiC31 
and Rv-BamHI-PhiC31 and cloned into the Pstl and 
BamHl sites of pBEVY-Nat giving the pBEVY-Nat- 
PhiC31. 



Table 1 Saccharomyces cerevisiae strains used in the current work 



S. cerevisiae strains 


Genotype 


Plasmid 


Source 


CEN.PK 102-5b 


Mat a, ura 3-52, his2A1, leu 2-3/1 12 




(van Dijken et al. 2000) 


CEN.PKc 


CEN.PK 102-5b 


[pYX012; pYX022; pYX042] 


This work 


TMB3400 


5. cerevisiae industrial strain 




(Almeida et al. 2009b) 


GSE16 


GS1 .1 1-26+ backcrossing with a segregant of Ethanol 
Red that is tolerant towards acetic acid; MATa/a 




(Demeke et al. 2013b) 


GSE16-YAP1 


GSE16-APJH:: TDH3p-YAP1-CYClt 




This work 


GSE16 - MCR1 


GSE16- YLR446W-l:.TPIp-MCR] 




This work 


GSE16 - MCR1-YAP1 


GSE16- YLR446W-1:.TPIp-MCR}, APJ1-V: TDH3p-YAPi-CYClt 




This work 


GSE16 - AAAPJ1 


GSE 1 6-APJ 1 /APJ 1 nattlVattL 




This work 
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Table 2 Plasmids used in the current work 



Plasmids Relevant features Origin 



n4?6GPF) 




IMi imhern et al 1QQS^ 


p426GPD-YAP1 


TDH3p-YAPl-CYCh 


This work 


pUG6 


kanMX flanked by loxP sites 


(Guldener et al. 1996) (EUROSCARF; 
accession number P301 14) 


YE-plac 112 KanR 


Multicopy, KanMX 


(Jeppsson et al. 2003) 


YE-plac 112 KanR-YAP1 


TDH3p-YAP1-CYCh 


This work 


pJET1,2 


Multicopy 


Thermo Scientific, Belgium 


pJET1,2-attB-KanMX-attP 


Multicopy, KanMX under TEFp 


This work 


p-miYAr 1 


I unjp-YAr i-lyl IX — rxaniviA atttvattr system, integrative 


This work 


pSTBIue-1 


multi-purpose cloning vector with dual kanamycin/ampicillin 
resistance. 


Novagen (EMD Millipore) 


pSTBIue-YLR446W 


pSTBIue-1 with the YLR446W gene cloned into the multiple 
cloning region 


This work 


pYX012 


Integrative; URA3; TPUp 


R&D System, Inc., Wiesbaden, D 


pYX012-LoxPkanMXLoxP 


pYX012, with the KanMX cassette flanked by loxP sites, deriving 

1 1 Ul 1 1 pUUO 


This work 


pYX01 2-LoxPKanMXLoxP-MCR1 


pYX012- LoxPKanMXLoxP with the MCR1 gene inserted in the 
MCS under the control of the TPI1 promoter 


This work 


pSTBIue-YLR446WA- LoxPKanMXLoxP-MCR1 


pSTBIue-YLR446W with the LoxPKanMXLoxP-MCR1 cassette, 
deriving from pYX012 LoxPKanMXLoxP-MCRI, inserted into 
the YLR446W sequence 


This work 


pSH65 


Centromeric plasmid, GAL1p-cre, ble r 


(Gueldener et al. 2002) (EUROSCARF; 
accession number P30122) 


pJETI ,2- attB-hph- attP 


Multicopy, hph under TEFp 


This work 


pBEVY-Nat-phiC31 


Multicopy, PhiC31 integrase 


This work 



Construction of S. cerevisiae strains 
S. cerevisiae GSE16-YAP1 

The open reading frame of the YAP1 gene was amplified 
from S. cerevisiae strain TMB3500 (Almeida et al. 2009b), 
using primers YAP1-F and YAP1-R. The amplicon was 
ligated into p426GPD (Mumberg et al. 1995) resulting in 
p426GPD-YAPl used for transformation of E. coli DH5a 
cells and followed by sequence verification. The cassette 
from p426GPD-YAPl was amplified using primers GPD- 
YAP1-F and CYClt-YAPl-R. The amplified cassette was 
ligated into pUG6 after restriction with Aatll and Pvull, 
resulting in pUG6-YAPl. Amplification of the two hom- 
ologous regions (HR) used for integration of the YAP1 
cassette into the APJ1 locus were performed from S. cere- 
visiae GSE16 (Demeke et al. 2013b) with primers HR1-F/ 
HR1-R and HR2-F/HR2-R. Amplicon for RH2 was first 
ligated into pUG6-YAPl and transformed into E. coli 
DH5a cells resulting in pYAPl-HR2. Amplicon for HR1 
was ligated into pYAPl-HR2 and transformed into E. coli 
DH5a cells resulting in pYAPl-HR2-HRl. The selection 
cassette (KanMX) flanked by the attB and attP sites was 
amplified from pJETl,2-attB-KanMX-attP using attBP-F 
and attBP-R. The amplicon was ligated into pYAPl-HR2- 
HR1 and transformed into E. coli NEB 5-alpha resulting in 
p-intYAPl. The nucleotide sequence of each amplicon 



was verified after every subsequent cloning step. S. cerevi- 
siae strain GSE16 was used for expression of the YAP1 
construct. In this strain, the expression cassette containing 
the transcription factor YAP1 was integrated by linea- 
rization of the integrative cassette using Aatll and Notl, 
followed by transformation using a DMSO-modified ver- 
sion of the LiAc method (Hill et al. 1991). The selection of 
colonies was done on YNB plates with 150 ug.mL~ genet- 
icin G418 (Sigma). Verification of the correct insertion 
{APJ1 locus) was done by sequencing using primers Ver. 
insl-F/Ver.insl-R and Ver.ins2-F/Ver.ins2-R. The result- 
ing strain was named GSE16-YAP1. 

S. cerevisiae GSE16-MCR1 

The open reading frame of MCR1 was amplified from 
CEN.PK 102-5b using primers MCR1-F and MCR1-R. 
The amplified DNA fragment (1544 bp) was cloned into 
pSTBlue (£coRV site) and used to transform E. coli DH5 
a cells. A 1.5 Kb £coRI fragment containing MCR1 was 
isolated from pSTBlue-MCRl and cloned into the MCS 
of the yeast integrative plasmid pYX012, resulting in 
pYX012-MCRl. In parallel, the ORF YLR446W was am- 
plified from S. cerevisiae CEN.PK 102-5b using primers 
YLR446W-F and YLR446W-R and cloned into pSTBlue 
(£coRV site) resulting in pSTBlue-YLR446W. The selection 
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cassette (KanMX) flanked by two LoxP sites was amplified 
from pUG6 using primers YLR446W Lox-F and Lox-R, 
and cloned into pYX012-MCRl plasmid (Kpril site). The 
expression cassette (LoxPKanMXLoxP + promoter + MCR1 
ORF + terminator) was amplified using primers YLR446W 
Lox-F and YLR446W TER-R and cloned into pSTBlue- 
YLR446W after restriction with Btgl and BsrGl (this 
double digestion allowed the removal of the inner part of 
the YLR446W gene (-601 bp)). Each amplicon was veri- 
fied by sequencing analysis after every subsequent cloning 
step. S. cerevisiae strain GSE16 was used for expression of 
the MCRl construct. The expression cassette containing 
MCR1 was integrated into GSE16 after PCR amplification 
using primers YLR446W-F and YLR446W-R. Transform- 
ation was carried out using a DMSO-modified version of 
the LiAc method (Hill et al. 1991). Correct insertion in the 
YLR446W locus was verified by PCR. The resulting strain 
was named GSE16-MCR1. The removal of the dominant 
marker (KanMX) was obtained by transforming GSE16- 
MCR1 with pSH65. 

S. cerevisiae GSE16-MCR1-YAP1 

The strain was constructed from GSE16-MCR1 by inte- 
gration of p-intYAPl previously digested with Aatll and 
Notl (see S. cerevisiae strain overexpressing YAP1). The 
correct integration in the APJ1 locus was verified by PCR. 

S. cerevisiae GSE16- AAAPJ1 

The selection cassettes (KanMX and hph) flanked by 
the attB and attP sites were amplified by PCR from 
pJETl,2-attB-KanMX-attP and pJETl,2- attB-hph- attP 
using primers Fw-APJl-Al and Rv-APJl-A2 with 50 bp 
homologues regions. Deletion of the two APJ1 alleles was 
carried out by integrating both selection cassettes into 
the APJ1 loci of GSE16 (Demeke et al. 2013b) using an 
adapted electroporation method (Benatuil et al. 2010). 
APJ1 double deletion colonies were selected from YPD 
plates with hygromycin (300 ug.mL 1 ) and geneticin 
(G418) (200 ug.rnL" 1 ), and checked by PCR using the 
primers Fw-APJl-check and Rv-APJl-check. To loop out 
the markers, the colonies were transformed with pBevy- 
Nat-phi31 and selected in YPD plates with 100 ug.mLT 1 of 
nourseotricin. These colonies were also checked on YPD 
hygromycin and YPD geneticin plates. The plasmids con- 
taining the integrase were lost by growing the colonies in 
YPD liquid medium overnight and transferred twice. The 
resulting strain was named GSE16-AAAPJ1. 

Spruce hydrolysate Spruce hydrolysate was obtained 
from SEKAB E-Technology AB (Ornskoldsvik, Sweden), 
and consisted of the non-detoxified liquid fraction of 
spruce after a pretreatment by SO2 catalyzed steam explo- 
sion. It is referred in this work as spruce hydrolysate and 
had the following sugar composition: 11 g.L -1 glucose, 



17 g.L 1 mannose, 4 g.L 1 galactose, and 10 g.L 1 xylose. 
All the batch fermentations were carried out with the 
same batch of spruce hydrolysate which was kept at 4°C. 

Anaerobic batch fermentations of spruce hydrolysate 

Inoculum was prepared by growing the cells overnight 
in 1 L shake flasks containing 100 mL defined mineral 
medium with vitamins (Verduyn et al. 1992) with glucose 
as carbon source (20 g.L -1 ) buffered with phthalate buffer 
(50 mM, pH 5.0). After centrifugation and a washing step 
with deionized water, the pellet was resuspended with 
20 mL of fermentation medium, and immediately used for 
inoculation of the fermenter. The fermentation medium 
consisted of 100% (v/v) spruce hydrolysate supplemented 
with 10 g.L 1 of xylose, 1 g.L 1 yeast extract, 0.5 g.L 1 
(NH 4 ) 2 HP0 4 and 0.025 g.L" 1 MgS0 4 . 7H 2 0. The pH of 
the hydrolysate was adjusted to 5.0 with 8 M KOH prior 
to supplementation. The fermentations were carried out 
in 1.2 L MultiFors fermenters with 0.5 L working volume. 
Temperature was maintained at 30°C and the pH was kept 
at 5.0 by addition of 3 M KOH. Oxygen free conditions 
were maintained by sparging N2 at 0.2 L.min~ , and the 
agitation was set to 600 rpm. Reactors were inoculated 
to a biomass concentration of 1 g.L" 1 (cdw). Cultivations 
were performed in biological duplicates for each inves- 
tigated strain. Specific conversion rates of furfural and 
HMF were calculated assuming a pseudo-steady state 
during the exponential growth on glucose. 

Anaerobic batch fermentations in mineral medium 

Inoculum was prepared by growing the cells overnight 
in 100 mL shake flasks containing 20 mL defined mineral 
medium with vitamins (Verduyn et al. 1992) with glucose 
as carbon source (20 g.L~ ) buffered with phthalate buffer 
(50 mM, pH 5.0). Fermentations were carried out in 
100 mL glass bottles with 40 mL working volume of the 
same mineral medium and vitamins supplemented with 
5 g.L~ of glucose and 20 g.L 1 of xylose as carbon and 
energy sources. Temperature was maintained at 30°C and 
stirring was set at ca. 160 rpm. Oxygen-limited conditions 
were obtained by sealing the bottles and sparging N 2 at 
0.2 L.min 1 for at least 5 minutes before inoculation. A 
cotton-filled syringe was inserted through the rubber stop- 
per using a needle to avoid accumulation of gas inside the 
bottles. Bottles were inoculated to a biomass concentra- 
tion of ca. 1 g.L~ (cdw). Cultivations were performed in 
biological duplicates for each investigated strain. 

Analysis of substrates and products Controlled vol- 
umes of samples were taken regularly for analysis. Biomass 
was followed by OD620 measurements during the length 
of the fermentations and determination of cell dry weight 
measurements were also done at time zero (just after 
inoculation) and at times 42 h and 92 h. For biomass 
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determination, the cell pellet from 5 mL culture was 
washed with distilled water and dried on Gelman filters 
(0 47 mm Supor-450, 0.45 urn) in a microwave oven 
(350 W) for 8 minutes. Ethanol, glycerol, acetic acid, HMF 
and furfural were analysed by high performance liquid 
chromatography (HPLC; Waters Corporation, MA, USA) 
using an Aminex HPX-87H column (Bio-Rad, CA, USA) 
at 65°C. The mobile phase was 5 mM sulphuric acid with 
a flow rate of 0.6 mL.min -1 . Analysis of glucose, mannose, 
xylose, galactose and xylitol was performed on a Shodex™ 
SP-0810 sugar column (Showa Denko K.K, Japan) at 85°C 
with water as mobile phase and 0.6 mL.mkT 1 flow rate. All 
compounds were detected with a refractive index detector 
(Shimadzu, Tokyo, Japan, and Waters 2414, MA, USA 
respectively). Yields were calculated based on HPLC 
measurements. 

Results 

Genetic engineering of the xylose consuming industrial 
strain GSE16 

It has previously been reported that, whereas overexpres- 
sion of the transcription factor YAP1 resulted in a strain 
with a faster sugar consumption rate when fermenting 
60% (v/v) of spruce hydrolysate, the same strain was se- 
verely inhibited at higher concentrations of the substrate 
(Alriksson et al. 2010). However the study was performed 
using laboratory strains with a sensitive genetic back- 
ground towards hydrolysate-derived inhibitors (Martin 
and Jonsson 2003). From these results and further studies 
carried out in our group (data not shown), the importance 
of the genetic background and the initial inhibitors con- 
centrations on the effect of YAP1 overexpression was 
further highlighted, and stressed the necessity of assessing 
the effect of gene modification under more industrial rele- 
vant conditions, that is using up-to-date engineered indus- 
trial strains and process-like fermentation conditions. 

In this study, overexpression cassettes containing the 
genes that encode the transcription factor YAP1 and the 
mitochondrial NADH-cytochrome b5 reductase MCRl 
were designed and integrated in the genome of the ro- 
bust industrial strain GSE16, either alone or in combin- 
ation, generating strains GSE16-YAP1, GSE16-MCR1 and 
GSE16-YAP1-MCR1 (Table 1). Simultaneous YAP1 inte- 
gration and APJl deletion were obtained by inserting the 
YAP1 overexpression cassette into the APJl locus. The 
MCRl cassette was targeted to the YLR446W locus, a 
gene whose deletion has not affected the yeast perform- 
ance during fermentation (Subtil and Boles 2012). A con- 
trol strain in which both alleles of APJl were deleted 
(GSE16-AAAPJ1) was included during the fermentations 
(Table 1). 

The effect of each genetic modification was assessed dur- 
ing anaerobic batch fermentations of undiluted spruce hy- 
drolysate. In order to have a medium composition similar 



to large scale lignocellulose-based fermentations, spruce 
hydrolysate with limited nutrient supplementation was 
used (See Materials and methods). The xylose con- 
centration was increased from 10 g.L~ to 20 g.L _1 to 
allow a longer period for analysis of the xylose consump- 
tion phase. The initial composition of the fermentation 
medium consisted therefore of 11 g.L~ glucose, 17 g.L~ 
mannose, 4 g.L 1 galactose, and 20 g.L 1 xylose. The follow- 
ing inhibitors could also be identified: acetate (3.7 g.L _1 ), 
HMF (0.96 g.L" 1 ) and furfural (0.78 g.L -1 ). The strain 
robustness was first assessed by comparing lag phase 
duration, specific growth rate, and specific HMF and 
furfural conversion rates. Next, we evaluated the fermen- 
tation performance of the different strains in terms of 
glucose and xylose consumption rates, ethanol production 
rate, and product distribution. The carbon dioxide profile 
of all strains during the batch fermentations is presented 
in Additional file 1: Figure SI. 

Impact of strain engineering on growth and furaldehyde 
conversion 

During the anaerobic batch fermentations of 100% spruce 
hydrolysate, and with an initial biomass of 1 g.L^cell dry 
weight (cdw), the control strain GSE16 showed a lag phase 
of 9 hours (Table 3). The strains overexpressing YAP1 
showed a consistent decrease in the duration of the lag 
phase, which lasted for 5.3 ± 0.6 hours. The lag phase dur- 
ation was not significantly affected by overexpression of 
MCRl while it was increased (by ~3 h) for the strain 
AAAPJ1 (p =0.01). The specific growth rate was also 
altered in all the strains with the exception of AAAPJ1. 
GSE16-YAP1 showed a specific growth rate of 0.21 h _1 , 
that is an increase of around 60% when compared to the 
control strain (p = 0.08). A similar increase was also 
displayed by GSE16-MCR1 (p = 0.15) (Table 3). No add- 
itional increase was observed by the strain combining 
both modifications. 

As for furaldehyde conversion, all strains exhibited a 
higher specific conversion rate (g.g cells -1 .!! -1 ) for furfural 
than for HMF (Figure 1). The control strain displayed a 



Table 3 Lag phase duration and maximum specific growth 
rate during batch fermentation of spruce hydrolysate 



Strain 




Lag phase 3 (h) 


Growth rate, \i b (h _1 ) 


GSE16 




8.96 ± 0.08 


0.1 3 ±0.02 


GSE16 - 


YAP1 


5.29 ±0.59 


0.21 ±0.03 


GSE16 - 


MCR1 


7.65 ± 1 .06 


0.20 ± 0.04 


GSE16 - 


YAP1- MCRl 


5.49 ± 0.83 


0.21 ±0.02 


GSE16 - 


AAAPJ1 


1 2.04 ± 0.48 


0.1 2 ±0.02 



Initial biomass concentration was 1 g.L -1 cdw. The values show the mean and 
standard deviation of two biological replicates. 

a The lag phase is defined as the time between inoculation and the onset of 
the increase on carbon dioxide production. b Growth rates were calculated 
from the carbon dioxide production rates during the exponential phase 
on glucose. 
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Figure 1 In vivo specific furfural (A) and HMF (B) conversion rates (g.g cells" h ) displayed by the different strains during fermentation of 
spruce hydrolysate. The values were calculated during the first two hours (furfural) and ten hours (HMF), considering the initial cdw (1 g.L~ ). The 
bars represent mean values of two biological duplicates, and the error bars indicate the standard deviation. 



furfural specific conversion rate of 0.16 g.g cells _1 .h~ . 
GSE16-YAP1 showed the highest furfural specific conver- 
sion rate among all the strains (two times higher than the 
control strain, p = 0.002). Conversion of furfural was also 
improved when MCR1 was overexpressed (Figure 1). Like- 
wise, in vivo conversion of HMF was enhanced by over- 
expression of YAP1 or MCR1. GSE16-YAP1 was able to 
convert the inhibitor around two times faster (0.08 g.g 
cells _1 .h _1 ) than the parental strain (p = 0.02), and a 
67% increase (p = 0.11) was observed for GSE16-MCR1 
(Figure 1). As observed for the growth rate, overexpres- 
sion of YAP1 in combination with MCR1 did not result in 
additional improvements neither in furfural nor in HMF 
conversion capacity. Neither was any significant difference 
in reduction capacity observed between the control strain 
and the APJ1 double deletion mutant. 

Impact of strain engineering on sugar metabolism 

Cell dry weight (cdw) measured after 48 and 92 h of 
fermentation, suggested that all strains had a similar 
increase in biomass (data not shown). However, accurate 
determination of biomass concentration in the presence 



of lignocellulosic hydrolysate was made difficult by the 
presence of solid particles in suspension in the medium or 
because of variations in the broth color (in the case of 
optical density measurement) over time. Therefore the 
strains were compared in terms of volumetric rates (of 
sugar and product formation) instead of specific rates. 

The volumetric consumption rate of glucose was around 
0.81 g.L _1 .h _1 for the strains overexpressing YAP1 or 
MCR1, which was at least 65% higher than the rate 
observed for the control strain (Table 4). And again, no 
additive effect was seen in the strains combining the over- 
expression of both genes. Deletion of APJ1 resulted in a 
marginal increase in glucose utilization rate. Our analyt- 
ical set-up did not allow accurate quantification of the 
mannose and galactose consumption rates. However, 
approximate determinations showed that mannose was 
consumed shortly after glucose was depleted from the 
medium. Co-consumption of mannose and xylose was 
observed during this phase. After mannose was exhausted 
from the medium, the relatively low amount of galactose 
(4 g.LT ) was co-consumed together with xylose (data not 
shown). In contrast to glucose, the xylose consumption 



Wallace-Salinas et al. AMB Express 2014, 4:56 
http://www.amb-express.eom/content/4/1/56 



Page 7 of 12 



Table 4 Volumetric consumption rate of glucose and xylose and volumetric production rate of ethanol during 
fermentation of spruce hydrolysate 



Strain 



Volumetric consumption and production rates (g.L \h ') 







Glucose 


Xylose 


Ethanol (glucose phase) 


Ethanol (xylose phase*) 


Ethanol (overall) 


GSE16 




0.49 


±0.01 


0.22 ±0.01 


0.60 ±0.1 3 


0.1 3 ±0.00 


0.27 ± 0.04 


GSE16- 


YAP1 


0.81 


±0.15 


0.1 3 ±0.00 


0.90 ± 0.00 


0.09 ± 0.02 


0.21 ±0.00 


GSE16- 


MCR1 


0.81 


±0.12 


0.24 ± 0.02 


0.84 ± 0.08 


0.13 ±0.00 


0.27 ±0.01 


GSE16- 


YAP1- MCR1 


0.84 


±0.12 


0.1 3 ±0.00 


0.90 ± 0.08 


0.08 ± 0.00 


0.20 ±0.01 


GSE16- 


AAAPJ1 


0.64 


±0.03 


0.1 9 ±0.02 


0.57 ± 0.05 


0.11 ±0.01 


0.28 ±0.01 



+ During this phase galactose is also consumed. 



rate was found to be negatively affected by YAP1 overex- 
pression. While the control strain consumed 0.22 g xylose. 
L _1 .h~ , a 40% reduction was observed for the two strains 
overexpressing the transcription factor (p < 0.03). Overex- 
pression of MCR1 alone and deletion of APJ1 did not have 
a significant effect on xylose utilization (Table 4). Since 
the deletion of both alleles of APJ1 in GSE16-AAAPJ1 did 
not affect xylose consumption, it is very unlikely that the 
effect observed in the K4W-overexpressing strains (carry- 
ing only one deleted allele) could be responsible for the 
lower xylose consumption rate observed in these strains. 

The volumetric ethanol production rate was calculated 
both during the glucose and the xylose consumption 
phases (Table 4). As expected from the differences between 
the rates of consumption of the two sugars, the ethanol 
production rate measured during the xylose phase was 
considerably lower than the one observed during the 
glucose phase for all strains. For the control strain, 
for example, the ethanol production rate on xylose de- 
creased by 78% as compared to the one measured during 
the glucose phase. When comparing the different strains 
on glucose, GSE16-YAP1 showed one of the fastest etha- 
nol production rates (0.90 g ethanol.L _1 .h _1 , p = 0.08). 
More generally, the ethanol production rate calculated for 
the new constructs was between 18% and 53% higher than 
the one obtained for the control strain, thereby matching 
the increase in the glucose consumption rate. An excep- 
tion to this correlation was GSE16- AAAPJ1, for which 
the ethanol rate during the glucose phase was the same as 
for the control strain (Table 4). On xylose, both YAP1- 
overexpressing strains exhibited a 40% reduction in the 
ethanol production rate when compared to GSE16. For 
GSE16-MCR1 and GSE16- AAAPJ1, the ethanol rate in 
the xylose phase was comparable with that of the control 
strain. When considering the total length of the fermen- 
tations (92 h), no substantial differences were observed 
between the strains in terms of the volumetric production 
rate of ethanol (p>0.1) (Table 4). 

All strains showed similar metabolite distributions, 
with the exception of xylitol which was around 2.4 times 
higher for the strains overexpressing YAP1 (p < 0.03) 
(Table 5). The ethanol yields at the end of the fermentation 



were very similar between all strains, and accounted for 
about 50-58% of the maximum theoretical yield (Table 5). 
These low ethanol yields are very likely a result of a sig- 
nificant rate of evaporation in the fermenter. As previously 
reported (Bengtsson et al. 2009; Wahlbom et al. 2001), the 
low boiling point of ethanol and the continuous sparging 
of the fermenters with nitrogen importantly affect the 
ethanol yields during prolonged fermentations. As men- 
tioned earlier, the biomass yield was also similar for all 
strains as was the glycerol yield. The acetate yield was very 
low for all the strains and less than 0.01 g/g in all the cases 
(data not shown). 

In summary, the deletion of APJ1 in GSE16 was not 
detrimental for the strain performance when fermenting 
spruce hydrolysate. While it was not possible to evaluate 
the positive benefit of this deletion with the set-up used 
for the fermentations (expected concentrations of ethanol 
lower than 3% v/v), APJ1 can be considered as a good 
candidate for integration of expression cassettes in strains 
to be used during very high gravity fermentations. In con- 
trast, YAP1 and MCR1 had a positive effect on glucose 
fermentation in undetoxified hydrolysate but the effect 
was not cumulative. Finally, YAP1 overexpression had an 
unexpected negative impact on xylose utilization. 

YAP1 and xylose consumption in mineral medium 

The positive effect of YAP1 in relation to resistance to in- 
hibitors was clearly shown by the improved fermentation 
rate of glucose during the anaerobic batch fermentations 
of spruce hydrolysate. The effect during the xylose phase, 
on the other hand, was not clear. In order to further 
investigate the effect of YAP1 overexpression on xylose 
consumption and to prevent any interference caused by 
the complex matrix of spruce hydrolysate, the xylose con- 
sumption rate of strain GSE16-YAP1 was evaluated during 
anaerobic fermentation of mineral medium containing 
5 g.L~ of glucose and 20 g.L~ xylose. Strain GSE16 
was included as a control. With an initial biomass of 
1 g.L~ (cdw), both strains consumed all the glucose 
within 5 hours. However, the negative effect of YAP1 over- 
expression on xylose utilization was also observed in 
mineral medium, indicating that such effect was not 
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Table 5 Ethanol, glycerol, biomass and xylitol yields per 
gram of consumed sugars in anaerobic batch 
fermentation of spruce hydrolysate 



Strain 



■ ethanol 



Yields 1 (g.g -1 ) 

^glycerol Ybj omass 



'xylitol 



GSE16 0.27 ±0.08 0.05 ±0.00 0.04 ±0.01 0.05 ±0.01 

GSE16 - YAP1 0.28±0.05 0.05 ±0.01 0.04±0.00 0.12±0.00 

GSE16- MCR1 0.26 ±0.00 0.06 ±0.00 0.04 ±0.00 0.05 ±0.01 

GSE16- YAP1-MCR1 0.25±0.00 0.06±0.00 0.04±0.00 0.12 ±0.01 

GSE16 - AAAPJ1 0.30 ±0.00 0.04 ±0.00 0.04 ±0.00 0.07 ±0.02 

Yields of ethanol, glycerol and biomass were calculated based on the total 
consumed sugars, while xylitol yields were calculated based on the consumed 
xylose only. The numbers reported are means ± standard deviation (n = 2). 

related to the presence of inhibitors. After 54 h, the con- 
trol strain GSE16 consumed 13.06 ± 0.46 g of xylose, while 
GSE16-YAP1 only consumed 3.17 ± 0.1 g xylose (Figure 2). 
As observed during the fermentations of hydrolysate, the 
final yields of biomass and extracellular metabolites were 
similar between the strains, except for xylitol. After 138 h, 
the xylitol yield of GSE16-YAP1 was 93% higher than that 
of the control strain (Table 6). 

Discussion 

In the present study two different genetic modifications 
that have been reported to improve strain tolerance in a 
laboratory strain background were introduced in the in- 
dustrial S. cerevisiae strain GSE16 and assessed during 



fermentation of undetoxified spruce hydrolysate. The ob- 
jective was to test whether the overexpression of the 
selected genes would still be relevant in a strain with a 
high robust genetic background, to identify the most 
promising gene among the two candidates, and to assess 
any putative additive or synergistic effects. 

During the glucose consumption phase, the overexpres- 
sion of YAP1 was shown to be relevant for inhibitor con- 
version, even in a robust genetic background. Similarly, 
the strain robustness in relation to spruce inhibitors was 
further increased when overexpressing MCR1, but no cu- 
mulative effect of combined overexpression was revealed. 
Slighdy increased inhibitor tolerance was observed when 
overexpressing YAP1 as compared to MCR1, which might 
be explained by the diversity of genes controlled by this 
transcription factor with relevant functions for detoxifica- 
tion (reductases, transporters, and oxidative stress-related 
enzymes among others) (Alriksson et al. 2010; Herrero 
et al. 2008; Toone and Jones 1999). Still, the explanation 
for the almost equivalent effect obtained by MCR1 overex- 
pression is less evident. Previous results showed that over- 
expression of MCR1 resulted in better growth in mineral 
medium supplemented with 12 g.L _1 acetic acid. (Signori 
et al., personal communication). Considering the role of 
the enzyme encoded by MCR1 in maintaining the antioxi- 
dant D-erythroascorbic acid (EASC) in its reduced form 
(Lee et al. 2001), it is plausible that higher concentrations 
of this antioxidant may help counteracting the oxidative 




30 

time (h) 

Figure 2 Profile of xylose consumption for GSE16 (squares) and GSE16-YAP1 (circles) in small vials with mineral medium and 5 g.L -1 
glucose + 20 g.L -1 xylose. Initial biomass was 1 g.L~' cdw. The experiment was carried out in biological duplicates. The figure shows the data of 
a representative profile for each strain with deviation <10%. 
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Table 6 Ethanol, glycerol, biomass, xylitol and acetate yields obtained during fermentations on mineral medium with 
5 g.L -1 glucose and 20 g.L -1 xylose 



Strain 






Yield (g.g -1 ) 








YEthanol 


^Glycerol 


V 

1 Biomass 


Yxylitol 


^Acetate 


GSE16 


045 ± 0.03 


0.04 ± 0.00 


0.09 ± 0.00 


0.03 ± 0.00 


0.01 ±0.00 


GSE16-YAP1 


0.45 ± 0.02 


0.04 ± 0.00 


0.10 ±0.00 


0.06 ± 0.00 


0.01 ±0.00 



Yields of ethanol, glycerol, biomass and acetate were calculated based on the total sugar consumed, while xylitol yields were calculated based on the consumed 
xylose only. The numbers reported are means ± standard deviation (n = 2). 



effect exerted by the acetic acid present in the hydrolysate 
(Semchyshyn et al. 2011). In fact, the effect seen by MCR1 
overexpression in the presence of acetic acid could be 
similar to the effect obtained with the biosynthesis of 
ascorbic acid (ASC) reported by Branduardi (Branduardi 
et al. 2007). In this study the authors showed that in 
yeasts, biosynthesis of ASC -a molecule very similar in 
structure and properties to D-erythroascorbic acid (EASC) 
- conferred increased resistance to H 2 0 2 , low pH and 
organic acids. Moreover, the higher in vivo conversion 
rates of furaldehydes displayed by GSE16-MCR1 also 
points towards an involvement of the enzyme in the 
reduction (by a yet unknown mechanism) of HMF and 
furfural probably into their less inhibitory alcohol forms 
(Liu et al. 2005). Given the specific location of this 
enzyme, this result supports previous observations that 
indicated a damaging effect of hydrolysate-derived inhibi- 
tors to the mitochondria (Allen et al. 2010; Nguyen et al. 
2014). To the best of our knowledge this is the first report 
of a positive effect of MCR1 overexpression on hydrolysate 
detoxification; and all together the results suggest that 
increased concentrations of NADH-cytochrome b5 reduc- 
tase can improve the resistance of yeast to hydrolysate 
inhibitors and therefore increase the ethanol production 
rate on glucose. 

On xylose, the much slower sugar consumption rate 
(as compared to glucose) in one of the best reported pen- 
tose fermenting industrial yeasts (Demeke et al. 2013b) 
and its derivatives emphasized that anaerobic pentose 
metabolism still requires further improvement, especially 
in the presence of lignocellulosic inhibitors. Recent com- 
parisons of metabolic profiles between the glucose and xy- 
lose consumption phases of S. cerevisiae with (Wang et al. 
2014; Ask et al. 2013) and without inhibitors (Bergdahl 
et al. 2012; Matsushika et al. 2013), highlighted the signifi- 
cant perturbations in the metabolic capacities of the yeast 
caused by xylose. In these metabolomics studies, such 
perturbations were ascribed to the depletion of key metab- 
olites in glycolysis and cofactors (among other import- 
ant metabolic variations) suggesting inefficient metabolic 
states such as carbon starvation, and diminished biosyn- 
thetic capacities (Bergdahl et al. 2012; Matsushika et al. 
2013). 

Our results also revealed a far more complex set of cel- 
lular responses deriving from the interactions between 



xylose metabolism and YAP1 overexpression. And al- 
though the evaluated physiological responses do not pro- 
vide enough information for explaining the decrease in 
xylose consumption in K4W-overexpressing strains, the 
integration of these observations with the results of previ- 
ous studies could give some hints for further analysis. 

First, higher concentrations of xylitol were obtained 
for ^Pi-overexpressing strains, which may result from 
the unspecific reduction of xylose by the different reduc- 
tases whose transcription is under the control of YAP1 
(Toone and Jones 1999) In fact, yAPi-overexpressing 
strains showed a higher xylose reductase activity than 
the control strain when cells extracts were used for reduc- 
tion of xylose (data not shown). Possible inhibition of xy- 
lose isomerase (XI) by xylitol was considered, although 
the used XI originates from C. phytofermentans, and has 
been shown to be much less inhibited by xylitol than other 
previously expressed XIs in S. cerevisiae (Brat et al. 2009). 
Moreover, the lower xylose consumption for GSE16-YAP1 
was observed from the beginning of the fermentation, i.e. 
when negligible levels of xylitol had been formed. Besides, 
the rate of consumption was almost constant during most 
of the process (except at the end of the fermentation when 
the low concentrations of xylose probably reduce the 
conversion rate for all strains). This implies that the 
consumption of xylose did not vary in relation to the in- 
crease of xylitol in the medium, i.e. xylitol was not inhibit- 
ing the conversion of xylose to xylulose by XI. 

A second aspect was connected to the "history" of 
GSE16. Demeke and co-workers (Demeke et al. 2013b) 
reported that the diverse and complex paths followed 
during the development of GSE16 included unknown 
mutations that appeared to be linked to the high capacity 
of the strain to ferment xylose, but that such mutations 
could also be correlated with a possible detrimental effect 
in terms of inhibitor tolerance (Demeke et al. 2013b). The 
authors suggested that this mutually exclusive phenotype 
(good xylose fermentation - bad inhibitor tolerance and 
vice versa) could be either causally or structurally linked, i. 
e. genes responsible for improved xylose utilization could 
be functionally connected with genes responsible for the 
slower growth; or in the other case, such genes may be 
located close to each other in the genome (Demeke et al. 
2013b). The results obtained in the current study with the 
YAP1 overexpressing strains also point towards a mutually 
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exclusive phenotype, but in this case a structural link does 
not seem likely since the deletion of the locus used for 
YAP1 integration (APJ1) did not cause a detrimental effect 
on xylose utilization. 

With the exception of xylitol, comparable values were 
obtained for biomass and extracellular metabolites yields 
between the different strains in hydrolysate fermenta- 
tions, which suggested that the overexpression of YAP1 
affected the rate of xylose metabolism but not the product 
distribution. As similar lower xylose uptake was obtained 
in mineral medium we could conclude that the negative 
interactions between xylose metabolism and overexpres- 
sion of the Yapl transcription factor were not dependent 
on the presence of inhibitors and their associated cell re- 
sponses. When considering available omics data about the 
effects of overexpressing YAP1 in S. cerevisiae in mineral 
medium, 17 transcripts were up-regulated by overexpres- 
sion of the transcription factor (DeRisi et al. 1997) while 
around 55 proteins were present in higher concentrations 
(Jun et al. 2012). Although the differences in cultivation 
conditions do not permit to make any conclusion towards 
particular genes of interest, these two studies suggest that 
the overexpression of YAP1 imposes a higher demand to 
the cells biosynthetic capacity. This, together with an 
impairment of biosynthetic capabilities on xylose (as 
presented in the metabolomics studies previously com- 
mented (Bergdahl et al. 2012; Matsushika et al. 2013)), and 
the high cell maintenance energy required during growth 
on this sugar (Feng and Zhao 2013) could explain the 
slower growth on xylose in YAP1 overexpressing strains. 

Nevertheless, the unexpected results seen with YAP1- 
overexpression during xylose assimilation require deeper 
analysis to further understand the biological responses 
that limit the development of robust xylose consum- 
ing strains. In this respect, we consider that the study 
at the molecular level of the cellular responses of 
YAP1 -overexpression in GSE16 and other xylose con- 
suming strains (for example expressing the fungal redox 
pathway) during glucose and xylose utilization would re- 
veal important insights about limiting steps for xylose 
metabolism. 

In conclusion, overexpression of the transcription 
factor YAP1 and the mitochondrial reductase MCR1 
in the already robust strain GSE16 resulted in an 
even faster hexose catabolism in the presence of spruce 
hydrolysate -derived inhibitors, but the effect was not 
cumulative. The improved phenotype of MCR1 overex- 
pression seems to be related, at least in part, to a faster 
furaldehyde reduction, indicating that this reductase may 
have a wider substrate range than previously reported. 
Unexpected reduced xylose fermentation rate was obser- 
ved in YAP1 overexpressing strains and further studies 
are needed to elucidate the mechanisms behind this 
observation. 
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